ON THE DYNAMICS OF A TIME-PERIODIC EQUATION 
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QIUDONG WANG 

Abstract. In this paper we use the second order equation 

d 2 q . 2^dq i 2 • 

— + (A - 79 J— -q + q = fJ,q smut 

as a demonstrative example to illustrate how to apply the analysis of [WO] and 
[WOkj to the studies of concrete equations. We prove, among many other things, 
that there are positive measure sets of parameters (A, 7, /x, u>) corresponding to the 
case of intersected (See Fig. 1(a)) and the case of separated (See Fig. 1(b)) stable 
and unstable manifold of the solution q(t) = 0, t € M respectively, so that the 
corresponding equations admit strange attractors with SRB measures. 



In the history of the theory of dynamical systems, ordinary differential equations 
have served as a source of inspirations and a test ground. There are mainly two ways 
in relating the studies of differential equations to the studies of maps, both originated 
from the work of H. Poincare. The first is to use return maps locally defined on 
Poincare sections for the studies of autonomous equations. The second is to use the 
globally defined time-T maps for the studies of time-periodic equations. 

The studies of periodically forced second order equations, such as van de Pol's 
equation, Dufling's equation and the equation for non-linear pendulums, have played 
substantial roles in shaping the chaos theory in modern times ( |Ar] . [D] , [Vj . [Levlj . 
[Le] . |GH] ). When a homoclinic solution is periodically perturbed, the stable and the 
unstable manifold of the perturbed saddle intersect each other within a certain range 
of forcing parameters, generating homoclinic tangles and chaotic dynamics (jP], [Sj, 
[M]). See Fig. 1(a). There are also forcing parameters for which the stable and the 
unstable manifold of the perturbed saddle are pulled apart. See Fig. 1(b). 





(a) (b) 
Fig. 1 The stable and the unstable manifold of a perturbed saddle. 
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With respect to the two scenarios of Fig. 1, there has been a (somewhat incorrect) 
perception that the scenario depicted in Fig. 1(a) is much more interesting than that 
in Fig. 1(b), mainly due to the fact that Fig. 1(a) induces complicated dynamics 
while Fig 1(b) appears simple. A computational procedure (the Melnikov's method) 
has been developed to distinguish the two scenarios. See Sects. 4.5 and 4.6 of [GHJ 
for an exposition on Melnikov's method and its applications to a list of periodically 
perturbed equations. Historically, homoclinic tangles of Fig. 1(a) have been a strong 
influence, yet until very recently in the analysis of concrete equations, little has been 
rigorously proved beyond the existence of embedded horseshoes. 

In |WUj and [WOkj . we have connected the studies of periodically perturbed equa- 
tions in the vicinity of a dissipative homoclinic solution to the studies of two specific 
class of maps. The main idea, motivated by |ASj . is to abandon the globally defined 
time-T maps, and in its stead to compute the return maps in the extended phase 
space around the unperturbed homoclinic loop. It has turned out that (1) the dy- 
namics of homoclinic tangles of Fig. 1(a) are those of a family of infinitely wrapped 
horseshoe maps [WOkj (see Sect. 11.21 for details), and (2) the dynamics of Fig. 1(b) 
are those of a family of rank one maps [WO] . With the return maps obtained in [WO] 
and in [WOkj, we are now able to apply many existing theories on maps to the studies 
of periodically perturbed equations with dissipative homoclinic saddle. Among the 
theories directly applied are the Newhouse theory [N], [PTj on homoclinic tangency; 
the theory of SRB measures [SI], [R], [Bo] ; the theory of Henon-like attractors [BC2J, 
[MV] , |BY] ; and the recent theory of rank one chaos [WYl] - [WY3j based on the 
theory of Benedicks and Carleson on strongly dissipative Henon maps [BC2] . This 
method has led us to many new results. 

In this paper we illustrate how to apply the analysis of [WOj and [WOkj to the 
studies of concrete equations. We use the equation 

o do 
(0.1) ^ + (\ - iq 2 )-± - q + q 3 = fiq 2 smut 

as a demonstrative example to prove, among many other things, that there exist 
positive measure sets of parameters (A, 7, uS) corresponding respectively to the dy- 
namics scenarios of Fig. 1(a) and Fig. 1(b), so that equation (10.11) admits strange 
attractors with SRB measures. We note that the unperturbed part of equation (10.11) 
has been studied previously in jHR] . 



I. Statement of results 
Let us re-write equation (10.11) as 

dq 

Tt =p 

'1-1) ^7- = —(A - iq 2 )p + q - q 3 + uq 2 sm9 

dt 

d6 

— — UJ. 
dt 



The space of (q,p, 6) is the extended phase space for equation (jO.ip where 6 G S 1 . 



1.1. Return maps in extended phase space. We call a saddle fixed point of 
a two dimensional autonomous system a dissipative saddle if the magnitude of the 
associated negative eigenvalue is larger than that of the positive eigenvalue. We call 
a saddle point a homoclinic saddle if it is approached by one solution, which we call a 
homoclinic solution, from both the positive and the negative directions of time. The 
autonomous part of equation (10 .ip is written as 

(1.2) |f + (W) |_ ? + ? 3. , 

where A > and 7 are parameters. Let (q,p) G M 2 be the phase space where p — 
We start with a previously existed result on (11.20 [HRj . 

Proposition 1.1 (Dissipative homoclinic saddle). There exists Xq > suffi- 
ciently small, such that for A G [0, Ao), there exists a 7a, |7a| < 10A such that for 
7 = 7a; 

(i) equation has a homoclinic solution for (q,p) = (0,0), which we denote as 

£x = {£x{t) = (a x (t)Mt)), 
satisfying a x (t) > for all t; 

(ii) for any given L > 0, there exists a K(L) independent of X, such that for all 

te[-L,L], ' 

\l x (t) - £ (t)\ < K(L)X. 

See also [WJ for an alternative proof of this proposition. Note that £o(t) is a 
homoclinic solution of equation (11.21) for A = 7 = 0. Let 



(1.3) a = -{V>* + Z+\), /3=-{VX r +Z-\). 

Then — a, (3 are the two eigenvalues of (0,0). Since —a + (3 = —X < 0, (0,0) is a 
dissipative saddle provided that A > 0. 

Let A > be sufficiently small. We regard the value of A as been fixed. We 
say that A G (0, Ao) is non-resonant ii a, (3 satisfy certain Diophantine non-resonance 
condition. Namely, 

(H) there exists di, di > so that for all n,m G Z + , 

\na — m(3\ > di(\n\ + \m\)~ d ' 2 . 

From this point on we fixed the value of A G (0, Ao) and assume that it satisfies 
(H). Let 7a be as in Proposition ll.il We introduce a new parameter p for 7 by letting 

(1.4) 7 = 7A + /i p. 
u>, p, p are now the parameters of equation (10.11) . 

We turn now to equation ( 11. ip . Let I = i\ be the homoclinic loop of Proposition 
II . II in the space of (q,p). We construct a small neighborhood of £ by taking the union 



of a small neighborhood U £ of (0, 0) and a small neighborhood D around i out of 
Ui E . See Fig. 2. Let a ± G U £ H D be the two line segments depicted in Fig. 2, both 
perpendicular to the homoclinic solution. In the space of (q,p,0) we denote 

U £ = U £ x S\ T> = DxS 1 

and let 




D 



Fig. 2 U £ , D and a ± . 

Let A/" : S + — > £~ be the maps induced by the solutions of (11.11) on U £ and 
M. : S _ — > S + be the maps induced by the solutions of ( 11. lft on D. See Fig. 3. 
Let '■= A/" o : S _ — > S _ be the return map defined by equation (11.11) . We 

compute JF = T w ^ p ^ following the steps of |WUj and [WOkj . 




Fig. 3 M and M. 

We now define the range of parameters for this paper. Let R w be sufficiently large, 
and -fi^, R p and i? M be such that 



Ru << R p << Rfj,. 
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We only consider parameters (u, p, p) inside of 

¥={(uj,p,p): ^(0,^), pG(0,i? p ), /16 (O,^ 1 )}. 

In the extended phase space (q,p,9), q = p = is a periodic orbit. We study 
equation (11. II) . denoting the stable and the unstable manifold of the solution (q,p) = 
(0, 0) as W s and W u respectively. First we have 

Theorem 1 (Surface of tangency). There exists a continuous function S*(uj,p) : 
(0, RJ) x (0, R- 1 ) -> (0, R p ) such that 

(a) W u nW s ^®if0<p< S*(w, n); 

(b) rnr = li/ s*(w, p)<p<R p . 

So under the surface S* defined by p = S*(u,p) is the scenario of homoclinic 
tangles of Fig. 1(a) and above this surface is the scenario of Fig. 1(b). S* is one of 
the surfaces (the lowest one of the three) depicted in Fig. 6 in Sect. 11.31 

1.2. Homoclinic tangles. Assume that (u,p,p) G P is beneath the surface S* of 
Theorem [U We focus on the set of solutions that stay inside of U £ U D for all time. 
This set of solutions is the homoclinic tangle in the vicinity of £\. If (uj,p,p) G P is 
beneath S*, then the return map T = J 7 ^^,^ is only partially defined on S~ because 
some of the images of £~ under M. would first hit on the wrong side of the local 
stable manifold of (q,p) = (0, 0) in S + , then sneak out of U £ . See Fig. 4. 




Fig. 4 Partial returns to X . 

We prove that T is an infinitely wrapped horseshoe map, the geometric structure 
of which is as follows. Take an annulus A = S 1 x I (This is S _ for J 7 ). We represent 
points in S 1 and / by using variables 9 and z respectively. We call the direction of 
9 the horizontal direction and the direction of z the vertical direction. To form an 
infinitely wrapped horseshoe map, which we denote as J 7 , we first divide A into two 
vertical strips, denoted as V and U. T : V — > A is defined on V but not on U. We 
compress V in the vertical direction and stretch it in the horizontal direction, making 
the image infinitely long towards both ends. Then we fold it and wrap it around the 
annulus A infinitely many times. See Fig. 5. 
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Fig. 5 Infinitely wrapped horseshoe maps. 

Let 

Of = {(M) e V : ^"(M) G V Vn > 0}, = n n > ^ n (^). 

The dynamic structure of the homolinic tangle in the vicinity of i\ is manifested 
in that of A^. Ajc- obviously contain a horseshoe of infinitely many symbols. This 
horseshoe covers Smale's horseshoe and all its variations. It is the one that resides 
inside all homoclinic tangles. 

The structure of Qj? and A^- depend sensitively on the location of the folded part of 
J-(V). If this part is deep inside of U, then the entire homoclinic tangle is reduced to 
one horseshoe of infinitely many symbols. If it is located inside of V, then the homo- 
clinic tangles are likely to have attracting periodic solutions, or sinks and observable 
chaos associated with non-degenerate transversal homoclinic tangency. In particular, 
we have 

Theorem 2 (Dynamics of homoclinic tangles). There exists an open domain in 
the plane of (u,p), which we denote as DH C (0, Rj) x (0,R p ), such that DH x 
(0, -R^ 1 ) is strictly under the surface S* in P; and for any fixed (u, p) G DH, the 
return maps T p = Tu, P ^, as a one parameter family in p, satisfy the follows: 

(i) There exists a sequence of p, accumulating at p = 0, which we denote as 

R- 1 > pP > P? > ■■■> ^ > P^ > ■■■>0 

such that for all p G \jpn,Pn \, T = : Ajr — > Ajr conjugates to a full shift of 
countably many symbols. 

(ii) Inside each of the intervals [/4h-i, Pn ], n > 0, there exists p so that T p admits 
stable periodic solutions. 

(Hi) There are also parameters in each of the intervals [p^+x, Pn ], n > 0, such that 
T p admits non- degenerate transversal homoclinic tangency of a dissipative saddle fixed 
point. 

Remarks: (1) All results stated so far in this subsection about homoclinic tangles 
are new. We are not aware of any previous results of the same kind in the analysis 
of any concrete time-periodic second order equations. 



(2) Theorem is much stronger than the mere existence of an embedded horse- 
shoe. For these parameters, the entire homoclinic tangle around l\ contains not only 
nothing less [S.j, but also nothing more, than a uniformly hyperbolic invariant set 
conjugating to a horseshoe of countably many symbols. Since the attractive basin of 
such horseshoe is of Lebesgue measure zero, these tangles are not observable if one 
adopt the probabilistic point of view that only events of positive Lebesgue measures 
are observable in phase space. 

(3) Sinks of Theorem [2](ii) are sinks of very strong contraction. They are not related 
to those obtained through Newhouse theory. The existence of Newhouse sinks follows 
from Theorem [2]^iii). 

(4) It also follows from combining Theorem [2](iii) with the result of |MVj that 
there are positive measure set of parameters p, such that the return maps T admit 
Henon-like attractors. In contrast to the horseshoes of Theorem [2]^i), these are chaos 
in observable forms (SRB measures) |MVj . [BY] . 

1.3. Invariant curves and rank one chaos. Assume that (u),p,p) G P is above 
the surface S* of Theorem HJ We are in the case of Fig. 1(b). The return maps T are 
now well-defined on S~, and they are not in anyway less interesting than those for 
homoclinic tangles of Section 11.21 in terms of chaos dynamics. For these parameters, 
not only the return maps are rank one maps, to which the theory of [WYlj and [WY2] 
directly apply, but also they fall naturally into the specific category of rank one maps 
previously studied by Wang and Young in [WY4] . For these parameters the dynamics 
of equation (II. ip around t\ are determined by the magnitude of the forcing frequency 
uj. When the forcing frequency uj is small, we obtain, for all p sufficiently small, an 
attracting tori in the extended phase space. In particular, we obtain an attracting tori 
consisting of quasi-periodic solutions for a set of p with positive Lebesgue density at 
/i = 0. As uj increases, the attracting tori is dis-integrated into isolated periodic sinks 
and saddles. Increasing the forcing frequency u further, the stable and the unstable 
manifold of these periodic saddles will fold, and intersect to create horseshoe and 
strange attractors. We have, in particular, that these are strange attractors with 
SRB measures. In particular, we have 

Theorem 3 (Invariant curves and quasi-periodic tori). There exists a po > 
and a continuous function Q : (p ,R p ) x (0,-R^ 1 ) — > (0, R^), such that for any given 
(u,p,p) G P that is on the left hand side of the surface Q defined by uj = Q(p,p), 
that is, these satisfying uj < Q(p,p), we have the follows: 

(1) The return map JF : £~ — > £~ admits a globally attracting, simply closed 
invariant curve; and the maps induced on this invariant curve is equivalent to a circle 
diffeomorphism. 

(2) There is an open set DI C (0, RJ) x (p Q ,R p ) in the (uj, p) -plane, such that 
DI x (0,-R" 1 ) is located on the left hand side of the surface Q. Let (uj,p) G DI be 
fixed and regard = Tu, P ,ii a one parameter family. Then there exists a set of p, 
of positive Lebesgue density at p = 0, such that the circle diffeomorphism induced on 
the attracting invariant curve above is with an irrational rotation number. 
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Remark: See Fig. 6 for Q inside of P. Similar results are previously obtained 
for periodically kicked limit cycles in |WY4j . |WY5j . |LWYj . and for periodically 
perturbed hyperbolic periodic solutions in |Lev2] . [Hj . 

We turn to the case of strange attractors and rank one chaos. 

Theorem 4 (Strange attractors and rank one chaos). Let S*(u,fi) be as in 
Theorem d and 

S(u,n) = (l + u*)S*(w,fi). 
Let (u,p,p) G P be such that 

S*(u,fi) < p < S(u,p). 

Then, 

(1) the return map T : S~ — > Yr admits a global attractor that is strange in the 
sense that it contains an embedded horseshoe provided that uj > 100; and 

(2) there exists an open domain DC C (100, Rj) x (0, R p ) in the (u, p) -plane, such 
that DC x (0, -R^" 1 ) is located in between S* and the surface S defined by p = S(u, p). 
Let (u>,p) G DC be fixed and regard = J~u>,p,n as a one parameter family. Then 
there exists a set of p of positive Lebesgue density at p = 0, such that admits a 
strange attractor with an ergodic SRB measure v. Furthermore, almost every point 
of S~ is generic with respect to v. 

Remarks: (1) Various dynamical scenarios of Theorems [10 are put together in 
Fig. 6. This figure is self-explanatory. Results similar to Theorem H] are previously 
obtained for periodically kicked limit cycles in [WY4j . |WY5j . [LWY] . and for certain 
slow-fast systems in |GWYj . 

p Q 



Attracting 
invarient . 
curve 



Transition to 
Chaos 




Fig. 6 Dynamical scenarios in parameter space. 

(2) Theorems [3] and H] remain valid if we change the forcing function pq 2 sin ut in 
O.ip to simply p sin ut. This case is covered by [WO] . 



(3) SRB measures are constructed initially for uniformly hyperbolic systems by 
Sinai, Ruelle and Bowen (|5ij. [Rj . [Bo] ) . They are observable objects representing 
statistical order in chaos. We refer to [Y] for a review on the theory of SRB measures. 
[WY1] contains not only the existence of SRB measures but also a comprehensive 
dynamical profile for the good maps, including geometric structures of the attractors 
and statistical properties. We have opted to limit our statement to SRB measures, 
but all aspects of that larger dynamical picture in fact apply. 

(4) Finally we note that the open domains DH, DI and DC in Theorems [20 are 
not in any sense small in size. One main restriction is that their respective product 
to (0,-R" 1 ) do not intersect Q, S and S*. This is essentially the only restriction for 
DI of Theorem [3J For DH and DC we also require that u is reasonably large. 



2. Derivation of the return maps 

In this section we derive the return maps T = NoM. of Sect. Il.ll for equation (11 .11) 
assuming Proposition ll.il In Sect. 12.11 we introduce coordinate changes to transform 
the equations into canonical forms, which we will use in Sect. 12.21 to compute the 
return maps. 



2.1. Equation in canonical forms. In Sect. I2.1I A we normalize the linear part. 
In Sect. 12. IB we linearize equation (11. ip on U £ , the e- neighborhood of the saddle 
point. In Sect. 12. 1C we introduce a set of new variables and derive the corresponding 
equations on D, a small neighborhood around the homoclinic solution l\ out oiU\ e . 

A. Normalizing the linear part Let A G (0, Ao) be fixed satisfying (H) in Sect. 
II. H and 7 = 7a + /ip where 7a is as in Proposition 11.11 We re- write equation (11.11) as 



dq 

Tt =P 
dp 

(2.1) — = -(A - 7a<? 2 )p + q - q 3 + n(pq 2 p + q 2 sin 9) 
d9 

— = UJ. 

dt 

To put the linear part of equation ( 12.11) in canonical form, we introduce new vari- 
ables (x, y) so that 

(2.2) q = x + ay, p = —ax + y, 

where a = |(A + VA 2 + 4) is as in (jl.3p in Sect. 11.11 In reverse we have 

(2.3) x= 1 (q-ap) y = r-j (aq + p). 

1 + cr 1 + or 
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The new equations for (x,y) are 
dx 

— = -ax + f(x, y) + n(A(x, y)p + C(x, y) sin 9) 
dy 

(2.4) — = f3y + g(x, y) + n(B(x, y)p + D(x, y) sin 0) 
at 

de 

— — u 
dt 

where (3 = a -1 is again as in (11.31) . and 

f( x , v) = - r 2 (ta(« + ay) 2 (y - ax) + (x + ay) 3 ) , 

1 + cr ' 

9{x, y) = z ■ 1 2 (ta(jc + ay) 2 (y - ax) + (x + ay) 3 ) ; 
1 + a 1 

A (x, y) = . a 2 (x + ayf(y - ax), B(x, y) = 1 (x + ayf(y - ax); 
1 + cr 1 + cr 

a —1 

c(^> y) = ^— — 2"( x + a ^) 2 ' jD ( a; ' ^) = — 2 ( x + a ^) 2 - 

1 + a z 1 + cr 

Observe that the functions /, g>, A, 5, C, D are all functions of x, y of order at least 
two. In comparison to the equations studied in |WOk] . equation f)2.4p is in a slightly 
different form. In the rest of this section we present the derivations of the return 
maps following mainly |WOk] . To avoid repetitive writings, we will refer the reader 
to the corresponding part of |WOk] for the details of a few technical proofs along the 
way. 

Notations (u,p,fi) are the three parameters of equation (11.11) . In what follows 
we replace /i by ln/i. /i £ (0,-R" 1 ) corresponds to ln/i £ (— oo, — lni?^). Denote 
p = (uj, p, In//), p £ P = (0, Ru) x (0, R p ) x (— oo, — lni? M ). e is not a parameter of 
equation f 1 1.1 1) but an auxiliary parameter representing the size of U e . In what follows 
e is a positive number sufficiently small satisfying 

(2.5) R^ » e- 1 » R p » Ru. 

The letter K is reserved as a generic constant, the value of which varies from line 
to line. K is allowed to depend on parameters u,p and e but not p,. We also use 
K(e) to make explicit when a constant is a dependent of e. Letter K standing alone 
represents a constant independent of both e and p,. 

The intended formula for the return maps would contain explicit terms and "error" 
terms, and we aim on C 3 -control on all error terms in this paper. To facilitate our 
presentation we adopt a specific conventions for indicating controls on magnitude. 
For a given constant, we write 0(1), 0(e) or O(p) to indicate that the magnitude of 
the constant is bounded by K, Ke or K(e)p>, respectively. For a function of a set V 
of variables and parameters on a specific domain, we write Oy(l),Oy(e) or CV(/-0 
to indicate that the C r -norm of the function on the specified domain is bounded by 
K, Ke or K(e)p,, respectively. We chose to specify the domain in the surrounding text 
rather than explicitly involving it in the notation. For example, Og. p (p) represents a 
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function of 9, the C 3 -nor nfl of which with respect to 9, to, p and In p is bounded above 
by K{e)p. 

B. Linearization on U £ . Let X, Y be such that 

x = X + P(X, F) + pP(X, Y, 9; p) 
y = Y + Q(X,Y) + pQ(X,Y,9;p) 

where P,Q,P,Q as functions of X and Y are real-analytic on \(X, Y)\ < 2e, and 
the values of these functions and their first derivatives with respect to X and Y 
at (X, Y) = (0,0) are all zero. As is explicitly indicated in ( 12.61) . P and Q are 
independent of 9 and p. We also assume that 

P(X, Y, 9 + 2tt; p) = P(X, Y, 9; p), Q(X, Y, 9 + 2n; p) = Q(X, Y, 9; p) 

are periodic of period 2n in 9 and they are also real-analytic with respect to 9 and p 
for all 9 E R and p E P. We have 

Proposition 2.1. There exists a small neighborhood U £ of (0,0) in (X,Y)-space, the 
size of which are completely determined by equation U.l\) and di,d 2 in (H), such that 
there exists an analytic coordinate transformation in the form of \2. 6}) on U e = U e x S 1 
that transforms equation {2.J$ into 



^ = - aX — = (3Y - = u 
dt ' dt ' dt 

Moreover, the C A -norms of P,Q, P,Q as functions of X,Y,9, p are all uniformly 
bounded from above by a constant K that is independent of both e and p on (X, Y) e 
U £ , 9 E E and p E P. 

Proof: This is a standard linearization result. For this proposition to hold we need 
the non-resonant assumption (H) and the fact that f,g,A,B,C,D are terms of order 
higher than one at (x, y) = (0,0). See for instance |CLSj for a proof. □ 

C. A canonical form around homoclinic loop. We derive a standard form 
for equation ( 12. 4p around the homoclinic loop of equation ( 11.21) outside of Ui £ . Let 
i = i\ be the homoclinic solution of Proposition IX.lt In (x, y)-space we write I as 

x = a(t), y = b(t) 

where t E M is the time. Let 



[u(t),v{t)) 



-i 



be the unit tangent vector of i at lit). We replace t by s to write this homoclinic 
loop as £(s) = (a(s),b(s)). Let 

e(s) = (v(s),-u(s)). 



We emphasize that ln^, not /x, is regarded as one of the forcing parameters. Derivatives with 
respect to /i and to In fi differ by a factor /i. 

11 



We now introduce new variables (s, z) through 

(x,y) = £{s) + ze(s). 

This is to say that 

(2.7) x = x(s, z) := a(s) + v(s)z, y = y(s, z) := b(s) — u(s)z. 

Let us further re-scale z by letting Z = {i~ l z. (s, Z, 9) are the variables we use in this 
paragraph. 

Before deriving the equations for (s, Z, 9), we first define the domain on which these 
equations are derived. Let — L~, L + be the times t = t\ hits Ui e . L ± are completely 
determined by i\ and e. Defined 

D = {(s,Z;p) : s G [-2L~ , 2L + ], \Z\ < Kx(e), p G P}. 

for some Ki(e) independent of /i. We have 

Proposition 2.2. The equation for (s, Z, 6) on D is in the form of 
dZ 

— = E{s)Z + (Hi(s)p + M 2 (s) sin^) + O a , z ,eM 
at 

(2-8) - = 1 + O s , z ,eM 

at 

d9 

— = U! 

dt 



where 
(2.9) 

and 
(2.10) 



E{s) = v 2 {s){-a + d x f{a{s), b(s))) + u 2 {s){f3 + d y g{a{s), b{s))) 
- u(s)v(s)(d y f(a(s), b(s)) + d x g(a(s), b(s))), 

Hi(s) = v{s)A(a(s), b(s)) - u{s)B(a(s), 6(s)); 
H 2 (s) = v(s)C(a(s), b(s)) - u(s)D(a(s), b(s)). 



Proof: The derivations of equation (I2.8P is identical to that of equation (2.13) in 
[WOkj . We note that Ki(s) for D is also defined explicitly in [WOkj . □ 

2.2. Derivation of Return Maps. In this subsection we derive the return map 
T = Af o M. by using Propositions 12.11 and 12.21 In Sect. I2.2I A we defined precisely 
the surfaces S 1 * 1 , and discuss the issue of coordinate conversions on these surfaces. In 
Sect. I2.2I B we compute Af by using Proposition 12.11 and in Sect. I2.2C we compute 
M. by using Proposition 12.21 T = Af o M. is then obtained in Sect. 12. 2D . 

A. Poincare sections E*. Let (X, Y, 9) be the phase variables of Proposition 12. II 
We define E 1 * 1 inside of U E fl D by letting 

E- = {(X,Y,0): Y = e, \X\ < 6 G S 1 }, 

and 

E+ = {(X,Y,9): X = e, \Y\ < Ki{e)n, 9 G S 1 }. 
Ki(e) is the same as in Proposition 12. 2[ 
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Let q G S + or £~. We represent q by using the (X, Y, #)-coordinates, for which we 
have X = e on S + and Y = e on S~. We can also use (s, Z, 0) of Proposition 12.21 to 
represent the same q. Before computing the return maps, we need to first attend two 
issues that are technical in nature. First, we need to derive the defining equations of 
S ± for (s, Z, 9). Second, we need to be able to change variables back and forth from 
(X,Y,6) to (s,Z,6) on E ± . Let 

(2.11) X = /i~ 1 X, Y = fT 1 Y. 

Proposition 2.3. (a) On S + , we have 

(i) s = —L~ + Ozep(^); and 

(ii) Y = (1 + 0{e))Z + O e , p (l) + OzfiM- 
(b) On we have 

(i) s = —L~ + Z fi, v {p); and 

(ii) Z = {1 + 0(e))X + O , P {1) + O^eM- 

Proof: The proof of (a) is identical to that of Lemma 3.1 and 3.3 of [WOkJ and the 
proof of (b) is identical to that of Lemma 3.4 of WOkJ. □ 

B. The induced map M : £+ -»• E~. For (X,Y,0) G S+ we have X = e/j" 1 by 
definition. Similarly, for (X, Y, 9) G X~ we have Y = eyT 1 . Denote a point on E + by 
using (Y, 9) and a point on X - by using (X, 9), and let 

(X X A) =JV(Y,0) 

for (Y,0) G £+. 

Proposition 2.4. We have for (Y,9) G S + ; 

Xi = (yU£ _1 )f 

(2.12) cj , 1, u 

Proof: Let T be the time it takes for the solution of the linearized equation of 
Proposition O from (e,Y,6) G £+ to get to (X l ,£,9 1 ) G £~. We have 

X l = ee- aT , e = Ye pT , 9 1 = 9 + uT, 

from which (I2.12p follows. □ 

C. The induced map M : £~ -> £+. Let H^s), z = 1, 2 be as in fl2TTU|) . In what 
follows, we write 

~L+ 



(2.13) 



We also write 



1 

H 2 (s) sin(# + us + wiT^o B(r)dr ds. 



(2.14) p L = e I- L - E ( s ) ds ^ p+ = e I E ( s ) ds _ 
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Note that for P L we integrate from s = —L to s = L + , while for P£ the integration 
starts from s = 0. First we have 

Lemma 2.1. 

P L no £0 « << 1, P+ ~ £ a >> 1. 

Proof: Same as Lemma 3.5 in [WOkj . □ 

For g = (s~,Z,9) G £~, the value of s~ is uniquely determined by that of 
(Z, 9) through Proposition 12.4( b) (i). So we can use (Z, 9) to represent q. Let 
(s(t), Z(t), 9(t)) be the solution of equation (12. 8p initiated at (s~,Z, 9), and t be 
the time (s(t), Z(i), 9{t)) hit £+. By definition M(q) = (s(t), Z(i),9(t)). In what 
follows we write 

Z = Z(t), 9 = 9(t). 

Proposition 2.5. Denote (Z,9) = M(Z,9). We have 

Z = Pt{pA L + Md)) + PlZ + O ZA M 
9 = 9 + u(L + + L-) + ZAp (fx). 

Proof: The same as that of Proposition 3.2 of |WOk] . □ 

D. The return map T = J\f o A4. We are now ready to compute the return map 
T = Af o M. : S _ — > We use (X, 9) to represent a point on E~ and denote 
(X 1 ,9 1 )=F(X,9). 

Proposition 2.6. The map T = J\f o M. : X~ — > £~ is given by 

X, = (^- 1 )f- 1 [(l + (9( £ ))P L + F(X^)]f 
(2 ' 16) 9 1= 9 + cu(L + + L-) + ^ In/i-yi + 0(e))P+ - ^ lnF(X, 9) + Xfi ,M 

where 

(2.17) 



F(X, 9) = (pA L + <j> L (9)) + PL{PiY\l + 0(e))X 
+ (Plr\l + Pl)O , p (1) + Ox,eM> 
and Pl,Pl an d 4>l{Q) are as i n Ii2.15\) and \2.1J$ . 
Proof: By using Proposition 12.51 and Proposition 12.3( b). we have 

Z = P L (1 + 0(e))X + P+{pA L + <j> L {9)) + P L O e , P {l) + Ox,eM 
9 = 9 + u(L + + L-) + O^M- 
Let Y be the Y-coordinate for (Z, 9), we have from Proposition 12.3( a). 

Y = (1 + £>(£))P+F(X, 9) 
where F(X, 9) is as in fl2~T7l) . We then obtain fl2~T6|) by using (I2TT2I) . □ 

14 



3. Proofs of Theorems CEH] 

Recall that i = i\ is the homoclinic solution of equation of Proposition 11.11 for 
equation (II. ip . which we write as £(s) = (a(s),b(s)). (u(s),v(s)) is the unit tangent 
vector of I at £{s). Also recall that Hj(s), i = 1, 2 are as in (I2.10p and E(s) is as in 
(TUD. Let 

/oo 
Il 1 (s)e-^ E( - T)dT ds 
-oo 

/oo 
U 2 (s)cos(ujs)e~fo E ^ dT ds 
-oo 



-oo 
»oo 



/oo 
U 2 (s) sm(ujs)e-ti E{T)dT ds. 
-oo 

First we need to estimate A, C(lo) and S{ui). 

Proposition 3.1. There exists a Ao(-R^), sufficiently small depending on R^, such 
that for all ui G (0, R u ) and all A G (0, Ao), we have 

(i) A > 1; and 

(ii) ^C 2 (u) + S 2 {u) > ( e -W + e^) _1 . 

Proposition 13.11 is proved in Section HI 

The return maps: To apply Proposition 12.61 we first fix a » 100/3. Then we 
pick a A G (0, Ao(-Rt,;)) satisfying the non-resonance assumption (H). We then chose 
Rp.e.R^ such that 

R^ » e' 1 » R p » A -1 . 
Recall that p = (u, p, In fi) G P represents the forcing parameter and 

£~ = {(9, X) : G E/(2ttZ), |X| < 1}. 
Let (^i,Xi) = T{Q, X) for (#, X) G £~ where T is from Proposition 12.61 We have 

6 1 = # + a- ^lnF(0,X,p) 

(3.2) P 

X! = b[F(0,X,p)]3 

where 

a = ^ ln/i" 1 + + L-) + ~ a ln(e(l + 

(3.3) P 
b = ( / xe- 1 )f- 1 [(l + 0( £ ))P+A i ]f 

and 

(3.4) F(0,X,p) = p + csin0 + kX + E(0,p) + C eiXiP (/i), 
in which 

(3.5) 



(A L )-\/ci + si 



k = P L (A L P+y\l + 0(s)) 
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and 

(3.6) E(9,p) = (A L P+)-\l + P L )O eiP (l). 

Note that in getting (13.21) we have changed 9 + ojL~ + c to 9 where c is such that 
tan c = C£ Sl. a, b, c, k and E(9,p) are as follows: 

(1) As parameters, b and a are functions of p = (u>, p, In p) and e. If we regard u, p 
and e as been fixed, and T p = T as a one parameter family, then b — » as \i — > 0. 
We can then think T as an unfolding of the ID maps 

f(0) = 9 + ti- ^ln(p + csinfl + E(fl,0)) 

P 

where a « — lo(5~ x In p. a — > +oo as p — >• 0. 

(2) a is a large number. But since it appears in the angular component we can 
module it by 2tt. Varying p, from R^ 1 > to zero is to run a over (ao, +oo) for some 
a ~ uP^lnR^. 

(3) c, k are independent of forcing parameters p but are functions of e. From 
Proposition 13.11 and Lemma 12.11 we have 

provided that e is sufficiently small. Recall that the value of e is selected after that 
of R w and R p . 

(4) From Lemma 12.11 and Proposition 13.11 we have 

E{9,p)^ £ ^~ 1 O 0iP (l). 
When e is sufficiently small, E(#, p) is a C r -small perturbation to p + c sin 9. 

Proof of Theorem [1} To prove theorem [TJ we fix the values of u,p and e and 

regard p G (0, R p ) as the parameter in T . Let 

F(9, p) = ¥(9, 0, p) = p + c sin 9 + E(9, p) + O e , p {p) 

and denote 

M(p) = min F(9,p). 
ees 1 

M(p) is a continuous function of p. 

We prove that M(p) is monotonic in p. For this it suffices to prove that for all 
fixed 9eS\ 

d p F(9,p) = 1 + d p E(9,p) + d p O e M > 1 - Kef - K(e)p > 0. 

Next we prove that IV s n W u = if and only if M(p) > 0. Let lf oc be the curve in 
£~ defined by X = 0. From (13.21) we know that y rn (l'i oc ) is well-defined for all n > 
if and only if M(p) > 0. We also know that W u n ^ u = if and only if F n (l? 0( ^ is 
well-defined for all n > because iy u came out of lf oc . 
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To finish the proof of Theorem [TJ it now suffices to observe that M(R p ) > and 
M(0) < because c ^ 0. The surface S* in P is defined by M(p) = 0. In fact we 
have 

(3.7) S*(u,fi) = c + 0(e + fi) 

from M(p) = 0. □ 

Proof of Theorem [2j The formulas (13.21) for the return maps are identical to 
these obtained in Sect. 4.1 of [WOk], and Proposition 13.11 implies that assumption 
(H2) in Sect. 4.1 of [WOkj hol dsE As sumption (HI) in Sect. 2.1 of [WOkj is (H). 
Consequently, Theorems 1-3 of [WOkj apply. A specific choice of the domain DH 
is defined in the paragraph entitled "specifications of parameters" of Sect. 4.1 of 
[WOkj . Theorem EKi) is Theorem 1 of jWOkj ; Theorem [2](ii) is Theorem 2 of [WOkj ; 
and Theorem [2](iii) is Theorem 3 of [WOkj . □ 

Proof of Theorem Let p > be such that 



p >2A-WC 2 (0) + S 2 (0) 

where C(0), 5(0) are the values of C{uj) and S(lj) at uj = 0. For (uj, p, p) G P satisfying 
P > Po, let 

M(u,p,fi) = minF(0,O,p). 

We have 

(3.8) M(0,p,p)>^p >0. 

5 

We now define Q by letting 

Q(p, p) = mm{u > : uj < 10' 5 M(uj, p, \i)cr\u)\ M(uj, p, p) > j^Po}- 

Q(p, p) is well defined because of ( 13.81) and the fact that M(uj,p, p) is continuous in 
uj. In the rest of this proof we assume (uj, p, p) G P is such that uj < Q(p, p). 

With (13.21) for T and the way Q(p, p) is defined above, our proof of Theorem [3] is 
the same as that of Theorem 1 of [WY4j . See Sect. 4.1 of jWY4j . For the existence of 
an attracting invariant curve, we observe that DT has a dominating splitting on 
(Sect. 4.1.1 of jWY4j ) . Theorem [3]( 2) then follows from Denjoy theory and a direct 
application of a Theorem of Herman [HeJ to the circle diffeomorphisms induced on 
the attracting invariant curves (Sect. 4.1.2 of [WY4j ) . 

To be more precise, we let v = (u, v) be a tangent vector of £~ at q and let 
s(v) = vu^ 1 for q = (9, X) 6 E". s(v) is the slope of v. Let Ch(q) be the collection of 
all v satisfying |s(v)| < and C v (q) be the collection of all v satisfying |s(v)| > 100. 
We first prove that 

(i) on E~, we have DT^C^q)) C C h (j T '(q)), and \DT q (v)\ > ||v|; 

(ii) on ^(S"), we have DF~ l {C v {q)) C C v {F{q)) and ID^ 1 ^)] > 100|v| for all 
v G C h (q) 



Though there is a difference in writing: in this paper we kept a factor p inside of F while in 
WOkj we gave it to b and a. With the current writing the proof of Theorem Q] is easier to present. 
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To prove (i) we compute DT by using (13.21) . Let (#i,Xi) = T(9, X), we have 
Mi Mi \ / i_ ig 



(3-9) ^= JL M = ^-iK^V^ 



96 



F i 



a 



where F = F(#, X,p) is as in (13.41) and 



_ = CCOsfl + ^"^(l) + O e ,X,a(/i) 

dW 



W/3 F^ 



< io- 4 



(i) then follows from b « p X/3 << 1 and 
(3.10) 

where (I3.10p is from lj < Q(p, p). To prove (ii) we use again (13 . 1 0[) and 

With both (i) and (ii), Theorem [3)^1) follows now from the standard graph trans- 
formation argument. Proof for Theorem [3)^2) is the same as the proof of Theorem 
l(a)(b) of |WY4j . See Sect. 4.1.2 of |WY4j . where T is the parameter a for here. □ 

Proof of Theorem |4j Theorem Hl(2) is proved by applying the theory of rank one 
maps of [WYT] and [WY2] to T of Q. The domain DC is such that 



202 y/C 2 (u) + S 2 {u) 396 y/C 2 {u) + S 2 {u 



<P< 



99 A r 101 A 

We also need to assume that u is sufficiently large for (w,p) G DC in order to prove 
that T^ is an admissible family. For more details see Section 6 of [WO] . 

Theorem H](l) is a specific case of Theorem 3 of [LW] . □ 

4. Proof of Proposition 13.11 

In this section we prove Proposition 13.11 We start with the conservative part of 
the autonomous equation ( 11.11) . that is, the equation 

(4.1) + ^ 

Denote p = -4. We write equation (14.1 p as 

( A ry\ d q dp 3 

(42) -dt=^ -dt =q - q - 

With a simple linear change of coordinates 

x = ^(q~p), y=^(q + p), 
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we write equation (14.11) in (x, y) as 

(4.3) 
Let 

a(t) 



doc 1 , ^ o 



— = y [x + y) . 

dt y 2 V yj 



St 



b{t) 



(l + e 2t)2> (1 + e 2 *) 2 ' 

(x, y) = (0,0) is a homoclinic saddle and (x,y) = (a(t),b(t)) is a homoclinic solution 
of equation (14.31) initiated at (\/2, v2). Denote 

{£(i) = ( a (t),&(t)), t G R}. 



It follows that 



(4.4) 



where 



u(t) 



-,2f 



3) 



VV - 3) 2 + (e- 2 * - 3)< 

e~ 2t - 3 
vV - 3) 2 + (e- 2t - 3) 1 



{u(t),v{t)) 

is the unit tangent vector of I at £(t). 
Let us re-write equation (14.31) as 



(4.5) 
where 

Let 

We have 
(4.6) 

Let 
(4.7) 



^ = -x + f{x,y) ^ = y + g(x,y) 



1 i 

f(x,y) = ^(x + y) 3 , g(x,y) = --{x + yf. 



E(t) = v 2 (t)(-l + d x f(a(t), &(*))) + u 2 (t)(l + d y g(a(t), &(*))) 
- u(t)v(t)(d v f(a(t), b(t)) + d x g{a{t), b(t))). 



E(t) 



3) s 



( e 2t _ 3)2 + ( e -2t _ 3)2 ^ + g 2t) 



12e 



2/ 



tf(s) = - / E(s)ds. 
Jo 

Lemma 4.1. For s G (— oo, oo), 

1 8e 2 *((l-3e 2s ) 2 + e 4g (e 2s -3) 2 ) 
(S) ~ 2 (e 2 * + 1)6 
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Proof: Let x = e 2t . We hav^l 

K(s) = 1 r Ul-3x) 2 -x 2 (x-3) 2 / _ 12s \ 
K > 2j l x{l-3x) 2 + x 2 {x-3) 2 V (l + x) 2 ) ' 

Observe that 

(1 - 3x) 2 + x 2 (x - 3) 2 = ((x - a) 2 + a 2 )((x - b) 2 + b 2 ) 
where a, b > satisfying a + 6 = 3, ab = \. We have 

1 r 2s (-x 4 + 6a; 3 -6a; + l)(x 2 -10a; + l) 
(S) ~ 2J 1 x{x + 1) 2 {{x - a) 2 + a 2 ){{x - b) 2 + b 2 ) 

e2s / 1 6 4x 3 - 18x 2 + 36a; - 6 



dx 



We then have 



2j x \x 1 + x ((x-a) 2 + a 2 )((x-b) 2 + b 2 ) 

1 _ _6_ + 2(x - a) + 2(s - 5) j ^ 
2 7 : \x 1 + x (x — a) 2 + a 2 (x — b) 2 + b 2 



1 8e 2s ((e 2s - a) 2 + a 2 )((e 2s - b) 2 + b 2 ) 



2 (e 2s + 1) 

1 8e 2s ((l-3e 2g ) 2 + e 4s (e 2g - 3) 2 ) 

2 n (e 2s + l) 6 

□ 

Through direct evaluation using Lemma 14.14 we obtain the values of the following 
definite integrals as 

/°° 16 
(u(s) + v(s))(b(s) + (a(s)) 2 (b(s) - a(s))e-^ E ^ dt ds = — . 

We have two more integrals to evaluate, and they are 

poo 

C{u) = / (u(s) +v(s))(b(s) + a(s)) 2 cos(ujs)e-f° mdt ds. 



DC 
DO 



(4.9) 

S(u) = I (u(s) + v(s))(b(s) + a(s)) 2 sin(^)e-^ i? W d ' c / s . 



Lemma 4.2. We have 



16V2n 2^2 + 



C(u) = — i J — , S(u 



2\ 



e -^um _|_ gjWTr 3 e -2 W7r -f e 2 L 

Proof: Using Lemma [4.11 we have 

C(uj) = 16V2I C , S(u) = 16y/2I s 



3 The author would like to thank Ali Oksasoglu for making this computation. 
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where 



r e 3 v l - e 2s ) 

I c = —rr, 7i- • cos (us) ■ as 

(e 2 * + l) 4 1 ; 



Let I — I c + il s . 

, (e 2s + l) 4 

We evaluate I as follows. On the z = x + iy plan, let 

t\ = {x, x : —00 — > 00}, £2 = {a^ + 7ri, x : 00 — > —00} 

and 

JlJ (e 2 * + l) 4 

We have 

/ f(z)dz = I, [ f(z)dz = e- u *I. 

So by the residue theorem, 

(1 + e~^)I = 2mRes(f(z)) 2 

Let t = z — y, we write f(z) as 

34/ 2* _j_ l) 



. 7T7. . 

2 



(1-e 



3 (5+u;i)i 1 p (3+wi)t 



16 t 4 (l + t + |t 2 + |t 3 + C(t 4 )) 4 ' 
We have 

e (5+Ht + e (3+<a,)t = i + ( 8 + 2^)t + (17 + 8^-cc> 2 )t 2 

+-(76 - 12w 2 + (51a; - cu 3 )i)t 3 + C(t 4 ). 
3 

We also have 

= l-4x + 10x 2 - 20x 3 + ■■■ 

(1+x) 4 

so 

(l+t + ^ + V = l-4t+—t 2 -8t 3 + 0(t 4 ). 
3 3 3 

Consequently, 

Res(f(z)) = ~e~^(8 - ±u(l + uj 2 )t), 



from which it follows that 



7 = 77- 8 - -w l + a; 2 i . 

8 1 + e- W7r v 3 v ; ; 
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In conclusion, we have 



U>(1 + u> ) 

Clio) = — j j — , o(^) = ■ — i 1 — • 

e -2 W7T + e2 W7r 3 e -2 OT + ez wn 



□ 



Proof of Proposition 13.11 We use subscript A to stress the role of A. We denote 
equation (12 .4p as (12.41) \. and the homoclinic solution £ x of Proposition 11.11 as 

x = a x (t), y = b x (t). 

Similarly, E x is defined through (I2.9p . and A\, C\{uj) and S\{uj) are defined through 
(13. ip . We need to estimate A x and C x {oj). 

In this section we have so far computed A x and C x , S x for A = 0. We have from 
and Lemma 

i a a 16 , N 16^71 

(4.10) A = — , C (lo] 



15 e ~2 U7T + ea 

We caution that a = /3 = 1 so for A = the saddle point (x, y) = (0, 0) is not 
dissipative. However, since —a + (3 — —A from (11.31) . (x,y) = (0,0) is a dissipative 
saddle for A > 0. We estimate A x and C\(uj) through A ,C (u)). 

Let U e be the e-neighborhood of (0, 0) in the (x, y)-plane, and s^(e) be the time £ x 
enters U £ . Similarly let — s^(e) be the time i x exits U £ . Then there exists Eq > and 
K independent of A such that for all A > sufficiently small, and all s G (s^(e ), +oo), 

(4.11) \u x {s)\ > 1-Ksq, \v x (s)\<Ke . 
Similarly we have, for all s G (— oo, — s^(eq)), 

(4.12) \v x (s)\ > 1 - Ke , \u x {s)\ < Ke . 

Both (14. lip and (14.121) follow directly from the fact that the local stable and the 
unstable manifold of equation (11.21) are tangent to the x and the y-axis respectively 
at (x,y) = (0,0), and their curvatures are bounded uniformly in A. From (14.111) and 
( I4.12p we have, for all s G (— oo, — s^(e )) U (s^(e ), +oo), 

(4-13) \E x (s)\ > \ 

provided that A are sufficiently small. Let 

A x , s -, s+ = f (u x (s) + v x (s))(b x (s) + a x (s)) 2 (b x (s) - a x (s))e-^ E ^ dt ds. 

J — s 

We have from (I4.13P that 

(4-14) \A X -A X>S -, S+ \<^L 

for all A that is sufficiently small provided that s~ > s^(eo), s + > s^(so). 
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Lets" = Sq(jE ), s + = Sq(jE ). From Proposition ll.l( ii). we have £\(— s~), i\(s + ) G 
U eo provided that A G (0, K~ 1 (eo)) where K(eo) > is a dependent of Sq. We also 
have £\(s) fl Ui £Q = for all s G (— s~, s + ), and this implies that for all s G (— s~, s + ) 

(4.15) \u x (s) - u (s)\, \v x {s) -v {s)\ < K(e )X. 
(I4.15P and Proposition ll.lf ii) together implies 

(4.16) \A x , s -,s+ - ^o,s-,s+ \ < K(e )X. 

That A\ > 1 follows now from P~T6lh (T4TT4j) and A > f|. 

Estimates on Ca(c^) is similar. Because C\{uS) is a function in u, so would be the 
constant K(e ) m (I4.16p . This is why A in this proposition is a dependent of R^. □ 
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